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Abstract 

The elliptic flow, V 2 , of muons from heavy-flavour hadron decays at forward rapidity (2.5 < y < 4) is 
measured in Pb-Pb collisions at = 2.76 TeV with the ALICE detector at the LHC. The scalar 
product, two- and four-particle Q cumulants and Lee-Yang zeros methods are used. The dependence 
of the V 2 of muons from heavy-flavour hadron decays on the collision centrality, in the range 0-40%, 
and on transverse momentum, pj, is studied in the interval 3 < px < 10 GeV/c. A positive V 2 is 
observed with the scalar product and two-particle Q cumulants in semi-central collisions (10-20% 
and 20^0% centrality classes) for the px interval from 3 to about 5 GeV/c. The V 2 magnitude tends 
to decrease towards more central collisions and with increasing px. It becomes compatible with zero 
in the interval 6 < px < 10 GeV/c. The results are compared to models describing the interaction of 
heavy quarks and open heavy-flavour hadrons with the high-density medium formed in high-energy 
heavy-ion collisions. The model calculations describe the measured V 2 within uncertainties. 
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1 Introduction 


Experiments with ultra-relativistic heavy-ion collisions aim at investigating the properties of strongly- 
interacting matter at very high temperatures and energy densities. Quantum Chromodynamics (QCD) 
calculations on the lattice predict, under these conditions, the formation of a Quark-Gluon Plasma (QGP), 
where colour confinement vanishes and chiral symmetry is partially restored |[T]-0|. Heavy quarks (charm 
and beauty) are created in initial hard-scattering processes on a time scale shorter than the QGP formation 
time. Subsequently, they interact with the medium constituents via inelastic |0 |7l and elastic lISUTOl 
processes. Therefore, heavy quarks are regarded as effective probes of the QGP properties. 

Heavy-quark energy loss due to in-medium interactions can be studied by means of the nuclear modifica¬ 
tion factor Raa, defined as the ratio of the yield of heavy-flavour particles measured in nucleus-nucleus 
(AA) collisions to that observed in proton-proton (pp) collisions scaled by the number of binary nucleon- 
nucleon collisions. The PHENIX and STAR Collaborations measured, in central Au-Au collisions at 
= 200 GeV, a strong suppression corresponding to a Raa of about 0.2-0.3 for heavy-flavour decay 
electrons at mid-rapidity (y) and transverse momentum pr > 5 GeV/c l[TTI - [r7l . A similar suppression 
was also measured by the STAR Collaboration for mid-rapidity D° mesons ifTSl . A significant suppres¬ 
sion was also observed by the PHENIX Collaboration at forward rapidity for muons from heavy-flavour 
hadron decays in central Cu-Cu collisions at = 200 GeV |[T^ . At the EHC, the AEICE Col¬ 
laboration reported a similar effect in central Pb-Pb collisions at y^NN = 2.76 TeV for D mesons at 
mid-rapidity Il20ll and muons from heavy-flavour hadron decays at forward rapidity |[^ in the interval 
2 < pt < 16 GeV/c and 4 < px < 10 GeV/c, respectively. The CMS Collaboration measured a signif¬ 
icant suppression of non-prompt J/t/r from beauty-hadron decays in the interval 6.5 < px < 30 GeV/c 
(3 < px < 30 GeV/c) and |y| < 2.4 (1.6 < |y| < 2.4) Il22ll23l . A first measurement of non-prompt J/t/r by 
the AEICE Collaboration at mid-rapidity (|y| < 0.8) and in the interval 4.5 < px < 10 GeV/c has been 
recently published ll^ . 


Eurther insights into the QGP evolution and the in-medium interactions can be gained from the study of 
the azimuthal anisotropy of particles carrying heavy quarks which, in contrast to light quarks, have expe¬ 
rienced the full system evolution. The study of azimuthal anisotropy is a field of intense experimental and 
theoretical investigations (see Il25]l and references therein). In non-central collisions, the initial spatial 
anisotropy of the overlap region, elongated in the direction perpendicular to the reaction plane, defined 
by the beam axis and the impact parameter of the collision, is converted into an anisotropy in momentum 
space through rescatterings ll2^ . Experimentally, the study of the particle azimuthal anisotropy is based 
on a Eourier expansion of azimuthal distributions given by: 


d^N 

dpjd(p 


1 dN 

2% dpx 


1 -h2 v„(px)cos[n((p -'P„)] 


n—\ 


( 1 ) 


where (p and px are the particle azimuthal angle and transverse momentum, respectively. The Eourier 
coefficients, v„, characterize the anisotropy of produced particles and is the azimuthal angle of the 
initial-state symmetry plane for the harmonic, introduced to account for the event-by-event fluctua¬ 
tions of the initial nucleon density profile. The second Eourier coefficient, V 2 , which can also be expressed 
as V 2 = (cos[2(<p —'^ 2 )]), is named elliptic flow. 

The V 2 of heavy-flavour hadrons is expected to provide information on the collective expansion and 
degree of thermalization of heavy quarks in the medium at low px (px <2 — 3 GeV/c). The participation 
of heavy quarks in the collective expansion is expected to give a positive V 2 ll26l. Moving towards 
intermediate px (3 < px < 6 GeV/c), the V 2 Eourier coefficient is also expected to be sensitive to the 
presence of recombination processes in the hadronization of heavy quarks ll27l l2^ . At high px (px > 
6 GeV/c), the V 2 measurement can constrain the path-length dependence of the in-medium parton energy 
loss, which becomes the dominant contribution to the azimuthal anisotropy and is also predicted to give 
a positive V 2 thus complementing the Raa measurement. 
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The PHENIX Collaboration reported a positive V 2 of heavy-flavour decay electrons at mid-rapidity in 
Au-Au collisions at -^^nn = 200 GeV, reaching a maximum value of about 0.15 at = 1.5 GeV/c in 
semi-central collisions imiBiiaii. A similar behaviour was also observed by the STAR Collabora¬ 
tion 132]]. Recently, a V 2 value significantly larger than zero was measured for D mesons at mid-rapidity 
in Pb-Pb collisions at = 2.76 TeV Il3^l34l . A complementary measurement at the same energy, 
provided by the heavy-flavour decay muon elliptic flow at forward rapidity (2.5 <y < 4), is of great 
interest in order to provide new constraints for models that implement the heavy-quark interactions with 
the medium. Einally, the measurement is also important for the interpretation of the J/t/r elliptic flow 
results at forward rapidity li^ in terms of a regeneration production from deconfined charm quarks in 
the medium. 

In this Eetter, we present the measurement of the elliptic flow of muons from heavy-flavour hadron decays 
at forward rapidity (2.5 < y < 4) in Pb-Pb collisions at .^^nn = 2.76 TeV recorded with the AEICE 
detector. The elliptic flow is measured using different methods: scalar product lIMIl . two- and four- 
particle Q cumulants |[37ll38l and Eee-Yang zeros |[39] - l4T]| . These methods exhibit different sensitivities 
to flow fluctuations and correlations not related to the azimuthal asymmetry in the initial geometry (non¬ 
flow effects). The V 2 coefficient is measured as a function of pj in the interval 3 < px < 10 GeV/c and in 
three centrality classes in the range 0-40%. The centrality dependence of V 2 is presented in the interval 
3 < pt < 10 GeV/c. 

The Eetter is organized as follows. The AEICE detector, with an emphasis on the muon spectrometer, 
and the data sample are presented in Section]^ The analysis details, the methods for the V 2 measurement, 
the inclusive muon V 2 determination, the procedure for the subtraction of the background of muons from 
decays of light-flavour hadrons and the study of systematic uncertainties, are described in Section 
The V 2 results for muons from heavy-flavour decays are presented in Section]^ and compared to model 
calculations in Section]^ Einally, conclusions are given in Section|^ 

2 ALICE experiment and data sample 

The AEICE detector is described in detail in Il42l l43l . The apparatus is composed of a set of central 
barrel detectors (pseudo-rapidity coverage |t]| < 0.9) located inside a solenoid magnet that generates a 
field of 0.5 T parallel to the beam direction, a muon spectrometer (—4 < rj < — 2.5^1 and a set of de¬ 
tectors for event characterization and triggering located in the forward and backward rj regions. The 
muon spectrometer consists of a passive front absorber made of carbon, concrete and steel, a beam 
shield, a 3 T-m dipole magnet, tracking chambers, a muon filter (iron wall) and trigger chambers. The 
muon tracking system is composed of five stations, each including two planes of cathod pad chambers, 
with the third station inside the dipole magnet. The muon tracking system is completed by four trig¬ 
ger planes of resistive plate chambers downstream of the iron wall, which absorbs hadrons that punch 
through the front absorber, as well as secondary particles produced inside it and low momentum muons 
(p < 4 GeV/c). Two scintillator arrays (VO) covering the pseudo-rapidity intervals —3.7 < T] < —1.7 
and 2.8 < T] < 5.1 are used for triggering, for collision centrality determination and for beam-induced 
background rejection. The Zero Degree Calorimeters (ZDC), located at 114 m from the centre of the de¬ 
tector on both sides, can detect spectator protons and neutrons and are also used for the offline rejection 
of beam-induced background and electromagnetic interactions. The Silicon Pixel Detector (SPD), that 
composes the two innermost layers of the Inner Tracking System (ITS), is used for the interaction vertex 
reconstruction. The Time Projection Chamber (TPC), which measures charged-particle tracks with full 
azimuthal coverage in |t]| < 0.9, is used in this analysis for the measurement of the reference particles 
(Section 3.1). 

Tn the ALICE reference frame, the muon spectrometer covers a negative rj range and consequently a negative y range. In 
the following, given that the colliding system is symmetric, the results are presented with a positive y notation. 
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The results presented in this Letter are obtained from the data sample recorded with ALICE during the 
2011 Pb-Pb run. The data were collected with a minimum-bias trigger requiring the coincidence of 
signals in the two VO arrays in synchronization with the passage of two crossing bunches. In addition, 
the recorded event sample was enriched with central and semi-central Pb-Pb collisions by applying 
thresholds, at the trigger level, on the VO signal amplitude. The beam-induced background (beam-gas 
interactions) was reduced by using the timing information from the VO and ZDC detectors. Lurthermore, 
a minimal energy deposit in the ZDC was required to reject the contribution from electromagnetic Pb-Pb 
interactions. Only events with a reconstructed primary vertex within ± 10 cm from the nominal position 
of the interaction vertex along the beam direction are analyzed. The Pb-Pb collisions are classified 
according to their degree of centrality by means of the sum of the amplitudes of the signals in the VO 
detector and the centrality classes are defined as percentiles of the total hadronic Pb-Pb cross section Il4^ . 
The analysis is carried out in three centrality classes: 0-10%, 10-20% and 20-40%. The analysed data 
sample corresponds to an integrated luminosity of 11.3 pb^^ in the 0-10% centrality class and of 3.5 pb^^ 
in the other two centrality classes. 


3 Data analysis 

The elliptic flow of muons from heavy-flavour hadron decays, 


of the inclusive muon elliptic flow, , 
pion and kaon decays (Sections 


is obtained from the measurement 
by subtracting the elliptic flow of muons from primary charged 
and: 


3.1 


3.4), as: 


u^HF ^2 
^2 = — 




( 2 ) 


where is the muon background fraction, defined as the ratio of the yield of muons from primary 

charged pion and kaon decays to that of inclusive muons. The measurement of the coefficient 

is carried out in the interval 3 < pj < 10 GeV/c in order to limit the systematic uncertainty on the 
subtraction of the muon background contribution. 


3.1 Track selection 

The selection criteria for particles of interest, muon tracks, are similar to those used in the previous 
analyses of pp collisions at y/s = 2.76 TeV and 7 TeV and Pb-Pb collisions at .^^nn = 2.76 TeV ll^l45]l . 
The tracks are required to be within the geometrical acceptance of the muon spectrometer, with — 4 < 
rj < —2.5 and 170° < Sabs < 178°, where Sabs is the polar angle measured at the end of the absorber. In 
order to improve the muon identification, a reconstructed track in the tracking chambers is required to 
match a track segment in the trigger chambers. This leads to a very efficient rejection of the background 
produced by charged hadrons, which are absorbed in the iron wall. Lurthermore, a cut on the product 
p ■ DC A of the track momentum p and distance of closest approach {DCA) to the primary vertex is applied 
to remove the beam-induced background tracks and fake tracks coming from the superposition of several 
particles crossing the muon spectrometer. This cut is set to 6a, where a is extracted from a Gaussian 
fit to the p DCA distribution measured in two intervals of 0abs. corresponding to different materials in 
the front absorber. After the cuts are applied, in the region pT > 3 GeV/c the residual background to 
heavy-flavour decay muons consists of muons from decays of primary charged pions and kaon^ and it 
amounts to 5-15%, depending on pT and on collision centrality (Section [T4] ). 

The mid-rapidity charged-particle tracks used to determine the flow vector Qn or the generating function 
(Section [T2] ) are called in the following reference particles. They are defined as tracks measured in the 
TPC in 1 77 1 < 0.8. These are required to have at least 70 associated space points out of the maximum of 

^Note that the contribution of muons from secondary light hadron decays produced inside the front absorber is negligible 
for pt > 3 GeV/c EH. 
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159, per degree of freedom (ndf) for the momentum fit in the range ;^^/ndf < 2 and a transverse 
momentum value in the interval 0.2 < pj < 5 GeV/c. Additionally, tracks are rejected if their distance 
of closest approach to the primary vertex is larger than 3 cm in the plane transverse to the beam direction 
or in the longitudinal direction. 


3.2 Flow analysis methods 

The elliptic flow measurement is carried out using various methods that have different sensitivities to 
flow fluctuations and non-flow effects Ii46l . Flow fluctuations are mainly due to event-by-event fluctu¬ 
ations of the initial density profile, while non-flow effects correspond to correlations not related to the 
azimuthal anisotropy in the initial state, such as resonance decays, jets and Bose-Einstein correlations 
between identical particles. It is worth mentioning that, in the present analysis, most of these non-flow 
effects are strongly suppressed by introducing an rj gap between reference particles and particles of in¬ 
terest HTl . In this analysis, the scalar product |[36l . two- and four-particle Q cumulants |[37l 1^ and 
Eee-Yang zeros Il39] - l4ll methods are employed. The description of these methods will be limited to 
the features specific to the present analysis. The following notations are introduced: 
refers to the measurement using the scalar product, and correspond to the 

ones using the two-particle Q cumulants and four-particle Q cumulants, while — Prod} 

and V 2 ^^^^^^{EYZ — Sum} are obtained using Eee-Yang zeros with product and sum generating func¬ 
tions. The superscripts p and p <— HE refer to inclusive muons and muons from heavy-flavour hadron 
decays, respectively. It is worth mentioning that these methods are more accurate than the standard event 
plane method, which yields a measurement lying between the event-averaged mean value and the root- 
mean-square value in the presence of flow fluctuations |[48ll49ll . Moreover, the multi-particle correlation 
methods (four-particle Q cumulants and Eee-Yang zeros) are less affected by non-flow correlations than 
two-particle correlation methods, but they cannot be used reliably when the muon flow magnitude is 
small and when the number of muons is small in the selected phase-space region e.g. in central and 
peripheral collisions, respectively llTTlI^ . Under these conditions, the scalar product and two-particle 
cumulant methods provide a V 2 measurement in a wider centrality range. 

The scalar product method llMll48l . derived from the standard event plane technique BSl . is based on 
the measurement of the flow vector Q„ computed from reference particles. In order to determine the 
elliptic flow, the Q 2 vector in a given event is expressed as: 


Qi 


N N 

7=1 7=1 


(3) 


where (pj is the particle azimuthal angle and N is the multiplicity of reference particles. 
With this method the 2"‘* harmonic coefficient is given by: 

V2{SP}= ^^2-Mr7^ ^ 

2V(ei-e1) 


(4) 


where the brackets in the numerator indicate the average over muons at forward rapidity, in all events. 
The vector Q 2 is calculated from Eq. ^ and the vector U 2 .i = (cos2<p,,sin2(p,) is the unit vector of the 
muon. In the denominator, each sample of reference particles used to compute Q 2 is divided into 
two sub-samples of same multiplicity in symmetrical rj intervals, —0.8 < rj < —0.5 and 0.5 < rj < 0.8, 
separated by a T] gap of one unit of pseudo-rapidity, labeled with the superscripts A and B and the 
brackets correspond to the average over events. 

The cumulant technique llT/l |38l is based on a cumulant expansion of multi-particle azimuthal correla¬ 
tions. Different order cumulants have different sensitivities to flow fluctuations. In the present analysis. 
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two- and four-particle cumulants are used to extract the muon elliptic flow. The results presented in 
the following are obtained from a direct calculation of multi-particle cumulants performed by using the 
2-cumulant technique lf38l . which is based on the moments of the magnitude of the flow vector Q 2 . It is 
worth mentioning that in this approach the cumulants are not biased by the interferences between various 
harmonics. The reference elliptic flow values V 2 evaluated from the order cumulant C 2 { 2 } and 4* 
order cumulant C 2 { 4 } with reference particles are given by V 2 { 2 } = y^C 2 { 2 } and V2{4} = C 2 { 4 }, 

respectively. Once the reference elliptic flow is estimated, the muon elliptic flow with respect to the 
reference elliptic flow is obtained from the 2 "‘^ and 4* order cumulants according to: 


V2{2} 


d2{2} 

1 / 2 { 2 } 


and V 2 { 4 } 


d2{4} 

E2{4}3’ 


(5) 


where <i 2 { 2 } and <i 2 { 4 } are the 2"‘^ and 4* order cumulants of selected muons Il38ll . 

The Lee-Yang zeros method |[39ti4TI relies on correlations involving all particles in the event. This is 
the limit of cumulants when the order of cumulants goes to infinity. The method is based on the location 
of the zeros in the complex plane, of a generating function of azimuthal correlations, which relates the 
position of the first minimum of the generating function to the magnitude of the reference elliptic flow 
V 2 defined as: 

l/2=(f;C0s[2((?)y-'P2)]) , (6) 

\ j—i ' events 

where M is the multiplicity of reference particles and the average is taken over all events. For this 
purpose, the following complex-valued generating function is evaluated as a function of a positive real 
variable r and few, typically five, equally spaced reference angles (LYZ — Prod method): 

=/n(l+«>cos[2((py-r?)])\ . (7) 

\y=l / events 

The first positive minimum of |G'^(/r)|, denoted as r^, allows one to estimate , which can be written 
as where 701 ~ 2.405 is the first root of the Bessel function. Once the first minimum 

is determined, the differential muon elliptic flow is estimated with respect to the reference flow 
detailed in flTI . Finally, the result is averaged over all t? angles. An alternative form of the generating 
function provided with the LYZ — Sum method is: 


G’^(/r) =/expr/r^cos[ 2 ((p 7 -r 7 )]'j\ . ( 8 ) 

The version of the method involving a product for the construction of the generation function (Eq. (|7]l) 
was designed to disentangle interferences between different harmonics, which is not the case with the 
generating function using a sum of the individual reference particle contributions. Both generating func¬ 
tions are used in this analysis. 

Note that, for all methods, autocorrelation effects are avoided because the particles (muons) used in the 
determination of the flow are not included in the estimation of the reference flow. 


3.3 Inclusive muon elliptic flow 

The elliptic flow of inclusive muons, is studied with two-particle correlation methods (scalar product 
and two-particle Q cumulants) in the centrality intervals 0-10%, 10-20% and 20-40%. In the 20-40% 
centrality interval, the multi-particle correlation methods (four-particle Q cumulants and Lee-Yang zeros) 
are also used. 
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Several sources of systematic uncertainty affecting the muon elliptic flow measurement are considered. 
These take into account the changes due to the variations of the reference particle selection criteria as 
in ll33|3l[5Ql, to allow us to check the robustness of the measurement. Since the collision impact 
parameter distribution could slightly depend on the observable used for the centrality determination, a 
systematic uncertainty is estimated by repeating the analysis using the number of clusters in the outer¬ 
most layer of the SPD and the number of tracks in the TPC as centrality estimators, instead of the VO 
signal amplitude. The systematic uncertainty due to the effect of TPC tracks from different Pb-Pb colli¬ 
sions piled-up in the same recorded event is estimated by applying a tighter cut to remove outliers in the 
multiplicity distribution of reference particles. This is done by requiring that the centrality values deter¬ 
mined using the VO signal amplitude and the number of TPC tracks do not differ by more than 5%. An 
additional systematic uncertainty specific to the scalar product is evaluated by varying the rj gap between 
the two sub-events from 1 to 0.8 T]-units (see Eq. Q and ll^ l. The various systematic uncertainties are 
added in quadrature. They tend to increase with increasing pj (see Eig.[^. A summary of the systematic 
uncertainties, in the interval 3 < px < 4.5 GeV/c, is presented in Table 1. 


V 2 analysis 

Source 

Systematic uncertainty (%) 

0-10% 

10-20% 

20^0% 

v^Tsp} 

Reference particles 

3 

1 

3 


Centrality selection 

6 

1 

4 


TPC pile-up 

2 

4 

2 


T] gap 

13 

1 

1 


Reference particles 

13 

3 

2 


Centrality selection 

14 

3 

6 


TPC pile-up 

8 

1 

4 

?{4) 

Reference particles 



10 


Centrality selection 



1 


TPC pile-up 



1 

V 2 {EYZ — Sum} 

Reference particles 



4 


Centrality selection 



7 


TPC pile-up 



2 

v^{EYZ-Prod} 

Reference particles 



2 


Centrality selection 



8 


TPC pile-up 



2 


Table 1: Systematic uncertainty sources affecting the inclusive muon elliptic flow measurement in the 0-10%, 
10-20% and 20^0% centrality classes for the interval 3 < px < 4.5 GeV/c. They are given as a percentage of the 
V 2 value. 

Eigureshows the px-differential muon elliptic flow (v^) in the 0-10%, 10-20% and 20^0% central¬ 
ity classes as obtained using the various methods. The values of slightly increase from central to 
semi-central collisions and this effect is more pronounced in the px interval 2 < px < 4 GeV/c. The 
two-particle correlation methods (scalar product and two-particle Q cumulants) give consistent results 
over the whole px range, indicating that these methods have a similar sensitivity to non-flow effect^ 
and in particular to flow fluctuations. A similar agreement is found when comparing the multi-particle 
correlation methods (four-particle Q cumulants and Eee-Yang zeros) to each other. No significant dif¬ 
ference between the results extracted with Eee-Yang zeros using either the sum or product generating 
function is seen, hence indicating that interferences between harmonics are negligible ||5T1 . The four- 
particle Q cumulants and Eee-Yang zeros are expected to be less affected by non-flow effects than scalar 

^Note that, in this analysis, most non-flow correlations are suppressed, even with two-particle correlation methods since 
reference particles and inclusive muons are separated by at least 1.7 p-units. However, it is worth mentioning that the main 
difference between the two methods is the p gap between the two sub-samples used to compute Q 2 (Eq. (|^) which also allows 
to partly remove non-flow effects. 
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Fig. 1: px-differential inclusive muon V 2 in 2.5 <y <4 and various centrality intervals, in Pb-Pb collisions at 
^snn = 2.76 TeV. The symbols are placed at the centre of the pj interval and, for visibility, the points from two- 
particle Q cumulants and Lee-Yang zeros with product generating function are shifted horizontally. The vertical 
error bars represent the statistical uncertainty, the horizontal error bars correspond to the width of the bin (not 
shown for the shifted data points) and the open boxes are the systematic uncertainties. The px intervals used 
with the Lee-Yang zeros method are different with respect to the other methods. Upper panels: results from two- 
particle correlation flow methods (scalar product and two-particle Q cumulants) in the 0-10% (left) and 10-20% 
(right) centrality intervals. Lower panels: results in the 20^0% centrality interval from two-particle correlation 
flow methods (scalar product and two-particle Q cumulants) and from four-particle Q cumulants (left), and from 
four-particle Q cumulants and Lee-Yang zeros (right). 


product or two-particle Q cumulants ||5^ . Moreover, since four-particle Q cumulants give comparable 
results as Lee-Yang zeros, one can conclude that non-flow correlations are almost completely removed 
at the order. Finally, the central values of obtained with four-particle Q cumulants or Lee-Yang 
zeros are systematically smaller than with two-particle correlation methods, although compatible within 
uncertainties. Such differences may indicate that initial fluctuations play a role in the development of the 
final momentum-space anisotropy. 

3.4 Muon background subtraction 

The subtraction of the muon background contribution to the measured requires an estimate of the 
elliptic flow of muons from charged pion and kaon decays, and of the background fraction, 

j-p-i^Tz.K i-ggg determination of the coefficient requires two steps. First, the pj- 

and T]-differential V 2 of charged particles measured in |t]| < 2.5 by the ATFAS Collaboration in Pb-Pb 
collisions 15^ and the px distributions of charged pions and kaons measured in |y| < 0.8 by the AFICE 
Collaboration in pp and Pb-Pb collisions ll5^ [55l are extrapolated to forward rapidity. Then, the px 
distributions of muons from charged pion and kaon decays, needed to estimate K are 
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generated according to a simulation taking into account the decay kinematics and the effect of the front 
absorber. 


The pt- and T]-differential elliptic flow of charged particles in |t]| < 2.5, v|*’, is extrapolated to forward 
rapidity using: 

V’?(PT,T?)=E(T])-vf(pT,2< |t]| <2.5), (9) 

where v|*’(pt,2 < |t]| < 2.5) is the measured charged-particle elliptic flow in 2 < |t]| < 2.5 with the 
event plane method. Since the V 2 *’(pt) measured by the ATLAS Collaboration is affected by statistical 
fluctuations, it is assumed that in the interval 10 < pj < 20 GeV/c, needed to simulate the decay muons up 
to pt = 10 GeV/c, v'^2 remains constant with a value given by the one measured in the interval 10 < pj < 
12 GeV/c. The extrapolation factor F(t]) is calculated by parameterizing the T]-differential v|*’ measured 
by the ATLAS Collaboration in various pj intervals with a second order polynomial. In the interval 
7 < Pt < 20 GeV/c, the ATLAS does not show a dependence on T] in Itj | < 2.5. Therefore, for px > 
7 GeV/c, F{r\) is computed as the average between a flat extrapolation function and the extrapolation 
factor obtained with the parabolic parameterization in 4 < px < 7 GeV/c. 


The mid-rapidity charged pion and kaon px distributions measured in Pb-Pb collisions are extrapolated 
to forward rapidity using the same strategy as in ll^ and summarized in the following. Assuming that 
the nuclear modification factor of charged pions and kaons in Pb-Pb collisions does not depend on 
rapidity up to y = 4 1121115^ . the px distributions of charged pions and kaons at forward rapidity can be 


expressed as: 


AN, 


n.K. 


PbPb 


dpxdy 


= (Taa) 


da, 


tc.K 

pp 


dpxdy 


[^AA (Pt)]>.=0, 


( 10 ) 


where (7aa) is the average nuclear overlap function in centrality classes under study, estimated as de¬ 
scribed in ||57]1 . The systematic uncertainty introduced by the assumption on will be discussed later. 
The rapidity extrapolation of the mid-rapidity pion and kaon px-differential cross sections measured in 
pp collisions 1121115^ is done according to: 


2_^,K 


d^a, 


pp 


dpxdy 


2^7r,K 


d^a, 


pp 


dpxdy 


J>>=0 




( 11 ) 


Gy being estimated from Monte-Carlo event generators (see 11211 for details). 


The elliptic flow of muons from charged pion and kaon decays, in 2.5 <y <4 and in various 

centrality classe^ is obtained by means of fast simulations using v|*’(77,px) given by Eq. (j^ and charged 
pion and kaon px distributions as obtained from Eq. The absorber effect is accounted for by 

rejecting the pions and kaons that do not decay within a distance corresponding to one interaction length 
from the beginning of the absorber. The simulation was repeated twice, considering that charged particles 
are either all pions or all kaons. The estimated decreases with increasing pj from about 0.085 at 

Px = 3 GeV/c to 0.035 at px = 10 GeV/c in the 10% most central collisions. In the 20-40% centrality 
interval, the values vary between 0.19 (px = 3 GeV/c) and 0.08 (px = 10 GeV/c) 


The background fraction, calculated as the ratio of the px-differential yield of muons from 

charged pion and kaon decays in 2.5 <y <4 obtained in the simulation to the measured px-differential 
yield of inclusive muons. It decreases as px increases, from about 12% (15%) at px = 3 GeV/c to 5% 
(7%) at Px = 10 GeV/c in the 0-10% (20-40%) centrality class. 

The systematic uncertainties affecting the estimated are summarized in Table 2. They originate 

from i) the method used to measure the charged-particle in ATLAS, ii) the T] and pj extrapolation of 


^The of muons from charged pion and kaon decays in the 20-40% centrality class is then obtained from the mean 

of the charged-particle V 2 in 20-30% and 30-40% centrality classes, with an additional systematic uncertainty provided by the 
difference with respect to the results in these two centrality classes. 
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V 2 '’ and iii) the treatment of the charged-particle in the fast simulation procedure. As the event plane 
method was used for the V 2 *’ measurement in ATLAS, the results range between the mean ((v^*’)) and 

R.M.S. {\J ((v|’’)^)) of the true v|*’ values due to fluctuations, depending on the event plane resolution 
which varies with the collision centrality 14^ . According to a Monte-Carlo Glauber model ll4^ . the ratio 
(V 2 )/(V 2 ) is expected to vary from about 1.06 to 1.15. Consequently, a conservative systematic uncer¬ 
tainty of 15% is applied to account for this bias and is propagated to The systematic uncertainty 

due to the T] extrapolation of is evaluated using several fit functions (first and third order polynomials, 
and Gaussian function) in the region pt <1 GeV/c, and for larger px values an additional systematic 
uncertainty due to the extrapolation procedure is considered. The latter is determined by comparing the 
results obtained with the two extrapolation functions used in the interval pj >1 GeV/c. The systematic 
uncertainty due to the assumption on in the region px > 10 GeV/c is estimated by varying be¬ 
tween 0 and the value in 10 < px < 12 GeV/c in the fast simulations. Such uncertainty affects mainly the 
high Px region (px > 7 GeV/c). Finally, the systematic uncertainty obtained by treating charged particles 
separately as pions and kaons is found to be negligible. The various systematic uncertainty sources are 
propagated to the estimated and added in quadrature. 


Source 

Systematic uncertainty (%) 

Input bias 

9 

V 2 ^ 77 extrapolation 

9-12 

high Px extrapolation 

13-15 

71 and K in fast simulations 

< 1 


Table 2: Systematic uncertainty sources affecting the estimated for the interval 3 < px < 10 GeV/c. They 

are stated as a percentage of the V 2 value. The given range reflects the dependence on the collision centrality. 

The systematic uncertainty on detailed in ||2T1 . includes the uncertainty on the generated px 

distributions of muons from charged pion and kaon decays, and the uncertainty on the measured inclu¬ 
sive muon Px distributions. The former originates from the input charged pion and kaon distributions, 
the rapidity extrapolation and the absorber effect. The systematic uncertainty on the measured inclusive 
muon yields contains the systematic uncertainty on detector response, residual mis-alignment and cen¬ 
trality dependence of the efficiency. This gives a total systematic uncertainty on about 21% 

in the interval 3 < px < 4.5 GeV/c with almost no dependence on the collision centrality. Finally, as 
done for the measurement of the heavy-flavour decay muon Raa 11211 . the systematic uncertainty due 
to the unknown suppression of charged particles at forward rapidity is calculated by varying 
from 0 to two times the estimated value. This corresponds to a variation of /?J'^(px) at forward rapidity 
from 0 up to two times [RJa ( 7 ^t)],v=o- This systematic uncertainty amounts to 10-30% in the interval 
3 < Px < 4.5 GeV/c, depending on the collision centrality and the flow analysis method. 

Finally, the systematic uncertainty on the elliptic flow of muons from heavy-flavour decays, V 2 con¬ 
tains two contributions: the systematic uncertainties on V 2 , V 2 and yM^^dc propagated according to 
the definition of given in Eq. (|^, and the systematic uncertainty due to the unknown suppression 

of charged particles at forward rapidity. The final systematic uncertainty on is obtained by adding 

in quadrature the two contributions. It amounts to about 12%-36% in the interval 3 < px < 4.5 GeV/c, 
depending on the collision centrality and the flow analysis method. 


4 Results 

Figure [^presents the px-differential elliptic flow of muons from heavy-flavour hadron decays, 
calculated with Eq. The results are shown for the 0-10% (upper, left), 10-20% (upper, right) and 
20-40% (bottom) centrality classes using the same flow methods as for the measurement of the inclusive 
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centrality intervals, in Pb-Pb collisions at ^snn = 2.76 TeV. The symbols are placed at the centre of the pj 
interval and, for visibility, the points from two-particle Q cumulants and Lee-Yang zeros with product generating 
function are shifted horizontally. The meaning of the symbols is the same as in Fig.[2 The horizontal error bars 
are not plotted for shifted data points. The pi intervals used with the Lee-Yang zeros method are different with 
respect to the other methods. Upper panels; results from two-particle correlation flow methods (scalar product 
and two-particle Q cumulants) in the 0-10% (left) and 10-20% (right) centrality intervals. Lower panels; results 
in the 20-40% centrality interval from two-particle correlation flow methods (scalar product and two-particle Q 
cumulants) and from four-particle Q cumulants (left), and from four-particle Q cumulants and Lee-Yang zeros 
(right). See the text for details. 


muon elliptic flow (Fig.^. When comparing the results to those obtained for inclusive muons (Fig.[T]), 
one can notice that and are similar due to the small background fraction (5% to 15%) in the 

Pj interval 3-10 GeV/c. The differences between the various methods are similar to those discussed for 
the measurement of the inclusive muon i.e. i) scalar product and two-particle Q cumulants give com¬ 
patible results, ii) consistent results are also found with four-particle Q cumulants and Fee-Yang zeros, 
and iii) the V 2 values extracted from these multi-particle correlation methods are smaller, although 
still compatible within uncertainties, than the ones obtained with two-particle correlation methods. As 
mentioned in Section 3.3 such differences are expected if initial-state fluctuations play a role in the 
development of the final momentum-space anisotropy. 


A positive is observed at intermediate px for the 20-40% and 10-20% centrality classes with 

a significance larger than 3a when combining statistical and systematic uncertainties. In the 20-40% 
centrality class, the values of the significance in the interval 3 < px < 4 GeV/c (4 < px < 5.5 GeV/c) 
are 4 a ( 3.2 a) and 4.3 a (3.8 a) with scalar product and two-particle Q cumulants, respectively. In the 
10-20% centrality class and in the interval 3 < px < 4.5 GeV/c, the values of the significance correspond 
to 4.4a both with scalar product and two-particle Q cumulants. This behaviour results from the interplay 
between the significant interaction of heavy quarks with the expanding medium and the path-length 
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Fig. 3: Elliptic flow of muons from heavy-flavour hadron decays as a function of the collision centrality in 2.5 < 
y < 4 and 3 < px < 10 GeV/c, in Pb-Pb collisions at = 2.76 TeV. The results are obtained with scalar 

product and two-particle Q cumulants. Vertical bars (open boxes) represent the statistical (systematic) uncertainty, 
the horizontal error bars correspond to the width of the centrality bin. For visibility, the points from scalar product 
are shifted horizontally and the horizontal error bars are not plotted 


dependence of in-medium parton energy loss ll29l[30l . The of muons from heavy-flavour hadron 

decays decreases with increasing px and becomes compatible with zero in the high pj region. 

Figurej^shows the centrality dependence of the pj-integrated (3 < px < 10 GeV/c) elliptic flow of muons 
from heavy-flavour hadron decays. It is investigated with scalar product and two-particle Q cumulants, 
which can be applied in a wider event-multiplicity (i.e. centrality) interval compared to multi-particle 
correlation methods. A significant decrease of the V 2 magnitude towards central collisions is observed. 
This is expected from the decrease of the initial spatial anisotropy from semi-central to central collisions. 

ALICE has measured the elliptic flow of prompt D mesons in |y| <0.8 in three centrality classes in the 
interval 0-50% with various two-particle correlation methods l|33j|34l. Similar trends as those reported 
here for muons from heavy-flavour decays are observed, although in different px and rapidity intervals. 
In particular, a positive V 2 was observed for D mesons in semi-central collisions in 2 < px < 6 GeV/c. 

The positive elliptic flow of muons from heavy-flavour hadron decays has been observed in a px interval 
from 3 to about 5 GeV/c where the charm contribution is expected to be dominant with respect to the 
beauty component according to perturbative QCD calculations ||2]- This measurement supports the 
interpretation of the J/tp positive V 2 at forward rapidity ll^ in terms of a significant contribution to J/tp 
production from recombination of flowing charm quarks in the deconfined medium. 

5 Comparison with models 

The results presented in this publication may constrain models describing the interactions of heavy 
quarks with the medium via elastic (collisional) and inelastic (radiative) processes, and in particular 
the parton energy loss dependence on the path-length within the medium. 
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Fig. 4: Left: px-differential elliptic flow of muons from heavy-flavour hadron decays in 2.5 < y < 4, in Pb-Pb 
collisions at = 2.76 TeV for the centrality class 20^0% compared to various transport model predictions: 
MC@sHQ H- EPOS ll59ll6Tll . TAMU igS and BAMPS ll63H65]l . The TAMU model is shown with a theoretical 
uncertainty band. Right: px-differential Raa of muons from heavy-flavour hadron decays for the centrality class 
0-10% from ED compared to the same models as for 

The elliptic flow coefficient and the nuclear modification factor of muons from heavy-flavour hadron de¬ 
cays ll^ are compared to the following three models. The MC@sHQ -i- EPOS transport model ll59ll treats 
the propagation of heavy quarks in the medium including collisional and radiative energy loss, within a 
3-1-1 dimensional fluid dynamical expansion based on the EPOS model if^lhTll . The hadronization of 
heavy quarks takes place at the transition temperature via recombination at low pj and fragmentation at 
intermediate and high pj. The final-state hadronic interactions are not included in the model. TAMU |[6^ 
is a transport model including only collisional processes via the Eangevin equation. The hydrodynamical 
expansion is constrained by pj spectra and elliptic flow data of light-flavour hadrons. The hadronization 
is modeled including a component of recombination of heavy quarks with light-flavour hadrons in the 
QGP The diffusion of heavy-flavour mesons in the hadronic phase is also included. BAMPS Il63l - l65l 
is a partonic transport model based on the Boltzmann approach to multi-parton scatterings. It includes 
collisional processes with a running strong coupling constant. The lack of radiative contributions is ac¬ 
counted for by scaling the binary cross section with a correction factor, tuned to describe the nuclear 
modification factor and elliptic flow results at RHIC energies. Vacuum fragmentation functions are used 
for the hadronization. 

Eigure (left) shows that the px-differential elliptic flow of muons from heavy-flavour hadron decays 
in the 20-40% centrality class is described reasonably well by the three models. However, the BAMPS 
model tends to slightly underestimate the Raa of muons from heavy-flavour decays in the 10% most cen¬ 
tral collisions, while the MC@sHQ-|-EPOS model tends to overestimate it. The TAMU model describes 
the Raa measurement over the entire pj interval within uncertainties and tends to slightly underestimates 
the V 2 measurement in the low pj region. This indicates that it is challenging to simultaneously de¬ 
scribe the strong suppression of high-px muons from heavy-flavour hadron decays in central collisions 
and the azimuthal anisotropy in semi-central collisions. Similar trends are also observed in the mid¬ 
rapidity region from the comparison of the Raa and V 2 of D mesons with model calculations ll34l . 

6 Conclusions 

In summary, we have reported on a measurement of the elliptic flow of muons from heavy-flavour hadron 
decays at forward rapidity in central and semi-central Pb-Pb collisions at .^^nn = 2.76 TeV with the 
AEICE detector at the EHC. 

Measurements have been carried out using several methods which exhibit different sensitivity to initial- 


13 

















Heavy-flavour decay muon V 2 in Pb-Pb collisions 


ALICE Collaboration 


state fluctuations and non-flow correlations. The systematic comparison of scalar product, two- and 
four-particle Q cumulants and Lee-Yang zeros helps in understanding the processes that build up the 
observed differences between two-particle correlation methods and multi-particle correlation methods 
and suggests that flow fluctuations are significant. 

The magnitude of the elliptic flow of muons from heavy-flavour hadron decays increases from central 
to semi-central collisions and decreases with increasing pj, becoming compatible with zero at high pj 
The results indicate a positive elliptic flow with the scalar product and two-particle Q cumulants in semi¬ 
central collisions (10-20% and 20-40% centrality classes) for the px interval from 3 to about 5 GeV/c 
with a significance larger than 3 a The data are described by transport model calculations within uncer¬ 
tainties, although a simultaneous description of /?aa and V 2 remains a challenge. The results reported in 
this Letter in various centrality classes may provide further important constraints to the models. 
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